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Abstract 

A method to derive homogeneous effective constitutive equations for periodically layered 
elastic media is proposed. The crucial and novel idea underlying the procedure is that the 
coefficients of the dynamic effective medium can be associated with the matrix logarithm of 
the propagator over a unit period. The effective homogeneous equations are shown to have 
the structure of a Willis material, characterized by anisotropic inertia and coupling between 
momentum and strain, in addition to effective elastic constants. Expressions are presented 
for the Willis material parameters which are formally valid at any frequency and horizontal 
wavenumber as long as the matrix logarithm is well defined. The general theory is exemplified 
for scalar SH motion. Low frequency, long wavelength expansions of the effective material 
parameters are also developed using a Magnus series and explicit estimates for the rate of 
convergence are derived. 

1 Introduction 

Elastic waves in periodically layered or continuous (functionally graded) elastic media of general 
anisotropy have been studied extensively by different methods. Among them is the sextic formalism 
of Stroh, which incorporates the elastodynamics equations into a first-order ordinary differential 
equation for the displacement-traction state vector with system matrix Q composed of material 
parameters [1]. The wave-field propagator matrix along the stratification direction y, M(y,0), is 
given by the Peano series of multiple integrals of products of Q (y) [2]. This is essentially a power 
series in distance-to-wavelength ratio, which is therefore particularly well-suited to tackling the 
problem of approximating a periodically stratified medium by an effective homogeneous medium. 
The Stroh formalism clarifies the meaning of zero-order homogenization, or static averaging, of a 
periodic medium by revealing the zero-order effective material parameters O 0] as nothing more 
than the matrix Q (y) integrated over the period T, which is the leading term of the logarithm 
of the Peano series for M(T, 0) [5]. Static averaging implies a non-dispersive effective medium. 
Generalization to a higher-order effective homogeneous medium, which must be dispersive, is less 
obvious. Its derivation is commonly based on the long-wave dispersion of the fundamental Bloch or 
Floquet branches. Their onset in arbitrary anisotropic periodically stratified media was analyzed in 
[5]; on this basis, the scalar- wave equation for a transversely isotropic dispersive effective medium 
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was modelled in [5l[6l[7] and in the subsequent literature (see e.g. [8] and its bibliography). A semi- 
analytical approach for general anisotropy [9] (see also [10]) fits the long- wave Floquet dispersion 
to statically averaged effective constants, such that the effective medium is seen as a "continuum 
of non-dispersive media" that are different for different frequency and propagation direction. 

In this paper a new method is proposed for finding dynamic effective constitutive equations at 
finite frequencies and wavelength. This is achieved by explicit construction of effective spatially 
constant material coefficients that reproduce exactly the monodromy matrix M(r, 0). The effec- 
tive constitutive theory is exact in the sense that it give the correct displacement-traction field at 
the unit-cell interfaces over arbitrary long distance of propagation. Two key steps distinguish 
the method advocated. First is the idea, based on the Floquet theorem, of defining the effective 
medium such that the sextic system of elastodynamics equations in this medium has the matrix of 
coefficients Qes equal to iK, where K is the Floquet wave number matrix with an exact definition 
zKT = lnM(T, 0). For the low-frequency long- wave range, a matrix logarithm lnM(r, 0) admits 
an expansion called the Magnus series [11]. Restricting it to the leading-order term leads to the 
statically averaged effective model (see above). Taking the next-order term(s) of the Magnus expan- 
sion for K reveals that, unless the variation of material properties over a period is symmetric, the 
above-defined dispersive effective medium cannot be fitted to the standard form of elastodynamic 
equations, in which the frequency dispersion and non-locality would be fully accounted by the de- 
pendence of effective density and elasticity on oj and kx- This motivates the second significant step 
in the present method which is identifying a class of constitutive models that does fit iK with the 
system matrix Qce of a homogeneous medium. We demonstrate by construction that the model 
described by the Willis constitutive relations with a dynamic stress-impulse coupling tensor [13^ [H] 
provides such a class of materials. Expansions of the Willis material coefficients based on the Mag- 
nus series for K are obtained in the low- frequency long- wave range where they are analytic in w, kx- 
Explicit estimates of the dependence on w, kx are found that enable closed- form asymptotics of the 
Willis coefficients with a desired accuracy. At the same time, the definition iKT = lnM(T, 0) and 
the fitting of the matrix Q^ff = iK to the effective Willis material is formally not restricted to the 
low-frequency long-wave range. 

A significant outcome of the proposed approach is that fairly explicit expressions are obtained 
for the effective material parameters. This is particularly evident for the example of SH (shear 
horizontal) waves discussed in detail in ^ In this regard the approach is distinct from that of |14j 
which leads to expressions for the parameters in the (spatial) transform domain. Note also that the 
results here apply to a single realization of the layered medium, no ensemble averaging is invoked. 
This point is discussed further in ^ 

The paper proceeds as follows. The problem is formulated in S}2] where the sextic formalism 
for periodically stratified media is outlined, the Floquet wave number matrix K introduced and its 
Magnus expansion examined (see also Appendix A) . The main results of the paper are derived in ^ 
It starts by observing that the matrix iK viewed as a sextic-system matrix Qcfj for a homogeneous 
effective medium cannot be associated with a medium from the class of anisotropic elastic materials 
but it does fit the Willis model. Using the ansatz that the sought effective medium is described by 
the spatially homogeneous equations for a Willis material, the corresponding system matrix Qcfr is 
constructed and equated with iK. Under certain assumptions, a prescription for unique definition 
of the Willis effective medium is put forward. The remainder of ^ discusses general properties 
of the Willis material parameters. The example of SH wave motion in a periodic structure is 
considered in ^ It illustrates the method for defining effective coefficients of the Willis model 
beyond the Magnus series expansion (which in its turn is detailed for SH waves in Appendix B). 
The obtained explicit expressions are used to solve a refiection-transmission problem at the interface 
of the effective medium. Conclusions are presented in ^ 
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2 Background 



2.1 Stroh formalism and the wave number matrix K 

We consider a Cartesian elastic medium with density p = p{y) and stiffness tensor cij^i = Cijki {y)- 
Basic notations used include tlie superscripts and * for transposition, Hermitian and complex 
conjugation, respectively, and T for the matrix with zero diagonal and identity off-diagonal blocks. 
Taking the Fourier transforms of the equilibrium and stress-strain equations 

= pUj, (Jij = CijklCkl, (1) 

in all variables except y leads to an ordinary differential problem for the quasi-plane modes with 
the phase factor ^K^^^-^'t) ^ where ui is the frequency and kx the wavenumber in an arbitrary chosen 
direction X orthogonal to Y (rotating X causes all 21 elastic constants Cijki to appear). Denote the 
unit vectors parallel to X and y by m and n so that x = m • r, y = n • r, and let A (y) and F (y) 
be the amplitudes of displacement u and traction na, respectively. Then Eqs. ([T]) combine into 
the Stroh system 

^r]{y) = Q{y)7]{y) (2) 
dy 

for the state vector r] incorporating A and F [1]. Taking it in the form tj = (A, iF)^ defines the 
6x6 system matrix as 

via the 3x3 blocks Nj of the Stroh matrix 

^ /Ni N2\ Ni = - (nn)-i (nm) = NJ, N2 = - (nn)"^ , 

' N3 = {mm) - {mn) (nn)"^ (nm) , ^ ' 

which is composed of the matrices with elements {nn)-^, = riiCijuni, {nm)-j^ = UiCijkimi = {mn)^,^ 

and [mm)-^ = miCijkimi (note that N2 is negative definite). The usual indicial symmetry of Cij^i 

used in (jH) leads to a Hamiltonian structure N = TN'^T and Q = TQ'^T of N and Q. Since cj, ky. 

and p, Cijki are real and hence Q (y) is imaginary, it follows that Q = — TQ^T. The latter identity 

on its own suffices to ensure energy conservation. Alternative definitions of r] and hence of Q may 

be chosen. In general, the eigenvalues of Q (y) are first-degree homogeneous functions of uj, kx, 

which implies absence of dispersion. 

Given the initial condition at some yo (= 0) , the solution to (l2|) is 77 (y) = M (y, 0) r] (0) , where 

M(y,0) is the 6x6 matricant evaluated by the Peano series [2] 

ry ry ryi ry ry^ cy^ 

M(y,0) = I+/ Q{yi)dyi+ Q{yi)dyJ Q (2/2) d2/2 + / / / ••. (5) 
Jo Jo Jo Jo Jo Jo 

Suppose now that p, Cijki and hence Q depend on y periodically with a period T. Denote 

y = y + nT, y = y modT G [0, T] , ?= | G [0, 1], i') = ^ '^^^^l'^^' 

where (...) is the static average over a unit cell. It is understood hereafter that the wave-path 
distance y includes a large enough number n of periods, which is when the present development is 
of interest. By ([5]), the matricant M (T, 0) over [0,T] , which is called the monodromy matrix, is 

M(r,o) = i + V" M^"^) =i + r(Q) + r2 /'Q(?)d?'rQ(?i)d?i + ... (7) 

^™=i Jo Jo 
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The wave number matrix K is introduced by denoting the monodromy matrix as 

M (T, 0) = exp (iKT) ^ iKT = In M (T, 0) . (8) 

In the fohowing, unless otherwise specified, K is understood as defined in the first Brihouin zone, 
which implies the zeroth Riemann sheet of Inz with a cut argz = itvr. Using ([8]), the matricant 
M (y, 0) can be written as 

M {y, 0) = M (y, 0) M (nT, 0) = L (y) exp {iKy) , (9) 

where L (y) = M (y, 0) exp (— iKy) with L (0) = L (T) = I. . Eq. ^2 represents the Floquet 
theorem. Denote the eigenvalues of M(r, 0) and K by e*^"^ and Ka (ce = 1--6), respectively. In 
the general case where M (T, 0) and K have six linear independent eigenvectors w^, the Floquet 
theorem implies that the wave field ry (y) = M (y, 0) rj (0) with the initial data expanded as rj (0) = 
Yla ^a^a takes the form 

^ (^) = ^ 1 <^"^" (y) ^he'^e (y) = L (y) ^ae'^'^y. (10) 

The identity Q = -TQ+T yields M'^ (y, 0) = TM+ (y, 0) T, which leads in turn to L^^ (y) = 
TL+ (y) T and 

K = TK+T = with K2,3 = K+3 (11) 

for K defined in the first Brillouin zone. If the unit-cell profile is symmetric, i.e. the variation of 
material properties within the period T is symmetric about the middle point so that Q (y) is even 
about y = T/2, then the above identities are complemented by M (T, 0) = TM^ (T, 0) T and hence 
K = TK'^T; so, with reference to ([H]), K is real. Thus 

K = TK^T = K* for a symmetric Q (y) . (12) 
2.2 Expansion of K in the Magnus series 

The logarithm of the monodromy matrix M (T, 0) can be expanded as a Magnus series |llj (see 
also Hg): 

^ 00 
iK = - InM (T, 0) = (Q) + ^ iK^"^^ with 

m=l 

zKW=Ti /'d?'r[Q(?),Q(?i)]d?i, (13) 
^ Jo Jo 

where [Q (x) , Q (y)] = Q (x) Q (y) — Q (y) Q (x) is a commutator of matrices depending on succes- 
sive integration variables. Each Magnus series term K^'"-' is a (m + l)-tuple integral of permutations 
of m nested commutators involving products of (m + 1) matrices Q (q) . A commutator-based form 
may be anticipated by noting that all K^'"-' for m > must vanish in the trivial case of a homo- 
geneous material with a constant Q = Qq and hence with M (T, 0) = exp (Qo^) . For practical 
calculations it is convenient to use the recursive formulas provided in [11] . In obvious contrast with 
the Peano expansion, the Magnus series converges in a limited range: the sufficient condition for 
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its convergence is (IIQII2) < t^/T, where || • is the matrix norm [T^. This condition imphes that 
the eigenvalues (uj'^,kx) of K, defined as continuous functions such that (0,0) = 0, satisfy 
the inequahty HRei^all <'k/T. 

The Magnus series for K is a low-frequency long-wave expansion. Actually uj and kx are two 
independent parameters for Q and hence for M and K. It is however essential that the dependence 
of 3x3 blocks of Q on a; and kx is homogeneous, see ([3]). Therefore the blocks of each mth term of 
the Peano and Magnus series are homogeneous polynomials of w, kx of degree one greater than the 
same block of the (m — l)th term. This is what enables introducing a single long-wave parameter 
e = kT with a suitably defined wavenumber A;, see (j55p . Thus the Magnus series (|13p is basically 
an expansion in powers of e. Taking small enough e enables its approximation by a finite number 
of terms. At the same time, it should be borne in mind that the Magnus series as an expansion of 
logarithm may converge relatively slowly. Explicit estimates expressed in terms of cj, kx and (N) 
which ensure a desired accuracy of truncating the Magnus series up to a given order are formulated 
in Appendix A. 

The structure of polynomial dependence of the Magnus series terms K^™) on kx and uj'^ is 

/ ki4Koo^4^ kx4^ \ 

-\kx{ki.^^+.-4') -ki4'^-u:-4'r 




(2)^/ kxikl^fUu:^^^) kl^f 

\kl{kl4^ + c^^af ) + a;^a(^) kx{kU?^ + M,^>')^ 

^(3) _ . ( kl{kl3.f^ + + cu^af kx{kl3.f + a;2af )) 

-\kx{kt.f+.-kl4Hu:^4') -kl{kl.f^ + .-.f^)-.^4'\ 



(14) 

etc., 



where the real matrices a^"^) in K^"^) are {m + l)-tuple integrals of appropriate commutators; for 
instance, a^^^ in K^^ are 

{a«, 4\ 4\ aj\ a(^)} = iT^'d?^'d?i|[Ni,Ni] + [N2,N3], [pl^n^], 

[Ni, N2] + [N2, , [Ns, Ni] + [N^ N3] , [Ni - nT pi] I . (15) 

The series terms K^"^) of odd and even order m are imaginary and real, respectively, and each term 
K(™-) on its own satisfies (jlip : therefore 

K^"^) = -K("*)* = -TK("*)'^T for odd m, K^"*) = K^™)* = TK^'^^^T for even m, (16) 

as taken into account in (jl4p . According to (jl2p . 

_ Q £qj, ■£ Q ^ -g gynimetric. (17) 

Note the pure dynamic imaginary terms proportional to ztiw™'^^, which appear in the diagonal 
blocks of the series terms K^™) of odd order m unless these are zero for a symmetric Q (y) by (jl7p . 

2.3 Dynamic homogenization 

According to the Floquet theorem ([9]), the wave field variation over a large distance y is characterized 
mainly by the function exp (iKy) (which is an exact wave field at y = riT). Formally exp (iKy) 
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with iKT = In M (T, 0) is a solution to Eq. ([2]) with the actual matrix of coefficients Q (y) replaced 
by the constant matrix iK. This motivates the concept of an effective homogeneous medium, whose 
material model admits the wave equation in the form ([2]) with a constant system matrix Qefr = iK. 
Confining the Magnus series (fT3]l to the zero-order term defines the statically averaged Q^.^ = (Q) 
which fits ([3]) with Neg = (N) and hence yields the non-dispersive effective density and stiffness in 
the well-known form p^'^^ = (p) and 

(nn)W = - (N2)-^ , {nmf^ = {N2)-' (Ni) , {mmf^ = (N3) - (nJ) (Ni) , (18) 

see [5]. It is evident that the statically averaged Q^g = (Q) completely ignores dynamic effects 
and is inadequate to describe waves at finite frequency over long propagation distance. Dynamic 
properties are realized by taking Qgff = beyond the zero-order term (Q) , see next Section. Note 
that, in contrast to (Q), a dispersive Qes = "iK generally depends on where the reference point 
y = of the period interval [0, T] is chosen. 

3 A dispersive effective medium with Qcff = iK 

3.1 The constitutive equations 

Our purpose is to take into account the full nature of the wave number matrix in Q^fj = iK. 
With this in mind, compare the structure of the dispersive effective matrix as given by the Magnus 
expansion Qcfj = (Q) + i X^^^i K^'") with that of Q (y) given by 1^. They differ in two ways. 
First, Qcff (7^ (Q)) is no longer imaginary, and hence the identity Qcff = — TQ^T which leads to 
(jlip (and ensures energy conservation) is no longer compatible with a Hamiltonian structure for 
Qeff. This is a well-known feature of dispersive models, see e.g. [T7]. The second, more significant, 
dissimilarity is that, by contrast to ([3]), Qcff has pure dynamic terms on the diagonal, already at 
the first-order zK^^^ see ()14p . This means that assuming dispersive density and elastic constants 
does not allow the constitutive relations of the dispersive effective medium to be written in the 
standard form of equations ([1]). Recalling that the upper rows of the sextic system ([2|) imply the 
traction-strain law, it is seen that the latter must be complemented by a purely dynamic term 
which implies different constitutive relations than those of the inhomogeneous medium itself. 

On this basis, following \13\ I14j. the equations of equilibrium and the constitutive relations ([T]) 
are replaced by the more general form proposed by Willis 

(Tij^i=Pj, (Jij = cf^i^ekl + SijrUr, Pq = SklqCkl + P^q^ ilr ■ (19) 

The vector p generalizes the normal notion of momentum density, and the elements of the Willis 
coupling tensor satisfy Sijk = Sjik by assumption, ensuring the symmetry of the stress tensor. A 
principal objective is to show that setting Qcff = iK leads inevitably to dispersive effective matrix 
density p'-^^^ and stiffness c\°^P and, on top of that, to the Willis form of the effective constitutive 
relations with stress-impulse coupling. 

3.2 The effective Wilhs medium 

Denote by Sn and Sm the matrices with components 

(Sn)jfc = n-iSijki (Sm)jfc = f^iSijk- (20) 
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The same derivation that led from ([T]) to the sextic system ([2|) with the coefficients ([3]) now leads 
from (I19p to ([2]) with the system matrix 



with L = S+nS'^^ + Nj^'^^+Sn + Sm + S+ = L+. The identity Qes = -TQ+gT is assumed 
in order to ensure that the effective medium is, like the inhomogeneous periodic medium, non- 
dissipative (energy conserving) . This implies hermiticity constraints on the complex- valued material 
parameters: c^^^^) = c^''^^+ , p^"^^ = p^''^^+ , S„ = -S*, = -S^,, where c^^*^) is the 6x6 
stiffness matrix in Voigt's notation. The coupling tensor Sijk is therefore purely imaginary. Note 

that the blocks N^''^^ of the effective Stroh matrix Neg = TN+jT TNjgT) consist of sub- 
matrices (nn)^^^^ , (nm)^'^^^^ , (mm)^^^^ built from c'^'^^-' according to the definition (j4]) but with 
(mnf''^ = (nm)(<'fif)+ so that Nf^ = Nf^+. 

Equating the matrix Qeff introduced in (I2ip to the matrix iK (w^, kx) with the block structure 
pT]) yields the blockwise equalities 

klNf^ - a;2(p(eff) _ S^Nf )S„) - ck^L = K3. (22) 

Identification of the effective parameters based on these identities is ambiguous given that they 
may depend upon both u and kx- A unique and arguably the simplest solution is obtained by first 
assuming a purely dynamic Sijk = Sijk (w). Then ([2^ ) and ^2% ) yield 

L^Sn = -K2 ^ (0) Ki (0) , (nn)('^) = -K2 \ A;^ (nm)^"^) = K^^Ki - K2 ^ (0) Ki (0) , (23) 

where the blocks Kj (^uj'^,kx) of K are taken at fe^, = as indicated by the notation Kj (0) = 
Kj (cj^,0) used here and subsequently. Note that the validity of ([23|) requires the block K2 to be 
invertible which is assumed in the following. For K defined by the Magnus series (jlSp . the conditions 
that ensure existence of Kg ^ in a certain low-frequency long- wave range and the estimates that 
enable truncation of its expansion are established in Appendix A. In view of (j23p . the remaining 
identity (j22]^ ) becomes 

kl (mm)(°fif) - cuV^'''^ - ujkx {S^ + S+ ) = K3 - K+K^^K^. (24) 

Pursuing the logical extension of the assumed dependence of the inertial coupling tensor S on 
frequency alone, we assume that the inertia tensor p^°^^ is also purely dynamic. This leads to a 
unique solution of (pi]) since p^'^^^ = p^°^^\lj) is found by setting k^ = 0, 

-a;(S™ + S+)= lim fc-ijKg - K+K^^Ki - (K3(0) - K+(0)K2 ^(O)Ki(O)) }. (25) 

The limit may be achieved in terms of derivatives of matrices Kj at k^ = 0, whose existence is 
guaranteed for instance within the range of convergence of the Magnus expansion. Accordingly, 
the solutions of (1241) are 



= K+(0)K2-H0)Ki(0) - K3(0), 
a;(S„, + S+ ) = -K's (0) - coS+K[ (0) - uKf (0) - uj^S+K'^ (0) S^, (26) 
kl {mmf^^ = K3 - K+K2 ^Ki + w2p(eff) + (S^ + S+ ) . 
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where K'j (0) = dKj {to'', A:,) /dk^^^^^. 

In summary, Eqs. (I23p and (I26p provide unique material properties for the Wihis model with 

c(°fif)(a;, A;^) = c(°fif)+, p(^*^)(a;) = p('=fif)+, S(l^) = -S*, (27) 

and Sijk = Sjik- The lack of dependence of the inertial parameters on kx means that non-local 
effects are confined to the elastic moduli c^''^^ The result reduces to non-dispersive statically 
averaged moduli of (|18p with p(°^) = (p)I and S = when fK is restricted to the zero-order term 

(Q) of (HID. 

Regarding computation of the Willis parameters from Eqs. (I23p and ()26p . it is assumed that the 
wavenumber matrix K (w^, k^) defined as iKT = InM (T, 0) is known either in the form of long- 
wave low-frequency series, or from the direct definition of matrix logarithm in some neighbourhood 
of a given point u, kx (see the example in fQ]). In particular, one may first evaluate the matricant 
M(T, 0) numerically and K then follows from the matrix logarithm. The matrix K'(0), required 
for the solution of Eq. (|26p 9. involves evaluating the derivative of lnM(T, 0) in kx- It may be 
expressed using either series or integral representation of In M as 



oo 1 

lnM = -V°° -(I-M)" = (M-I) / [x(m-l)+l]-^dx, 

^n=l n Jq 



which leads to 

n-l 



iK'{0)T = -(-1)""^ Yl A^M'(0)A"-^-J' = [ {xA + I]-^M'(0)[a;A + I]-^dx, 
n=l j=o •'^ 

with A = M (T, —\ ai kx = 0, and the derivative of the matricant itself is [2] 

aM(r, 0) '■^ 



(28) 



(29) 



M'(0) = 



dkx 



(30) 



The sufficient conditions for the range of validity of the above series and integral definitions of 
K'(0) are specified in Appendix A. 



3.3 Discussion 

3.3.1 The Willis equation and its inertial quantities 

Consider the above results (j23p and (j26p in more detail. Anisotropic density and the coupling 
coefficients that relate particle momentum and stress are unknown in "standard" models of solids. 
Here they appear as inevitable ingredients of a model that replaces periodic spatial inhomogeneity 
with a spatially homogeneous but dispersive and non-local theory. The departure from normal 
elasticity is evident from the equations of motion for the displacement that follows from (jl9p . 

4^u'^i,ik + {Siji - Siij)ui^i - pfpili = 0. (31) 

This in turn leads to an energy conservation equation of the form U +div f = where the real- valued 
energy density and flux vector are 

U = \4'jkhi,kuli + ^pfi^^iiiu*, fi = -^{Siji - Siij)uiu* - Re(c5J;^u/,fcU*) . (32) 
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In order to gain some insight into these new dynamic terms, consider a layered transversely 
isotropic medium, with the principal axis along n identified as the 2-direction. The effective density 
is of the form p^'^^^ = diatg{p^yl^\ p'^2\ Pu^^) the only non-zero elements of the coupling tensor 
(up to symmetries Sijk = Sjik) are 5'ii2 = 5'332, 5*211 = 5'233 and 5222- Only one combination of 
the three independent coupling elements has impact on the equations of motion, 

(^fjkl^l,ik + ('S'112 - 5'21l)u2J - Pll^Uj =0, j = 1, 3, 

(33) 

(^f2kl'^l,ik + (5211 - 5112) (mi,1 + 1X3,3) - P22^U2 = 0, 

and it has no influence on pure SH wave motion (polarized in the plane orthogonal to n). Long- 
wave expansions of p~^^^ and Sijk are presented in (pl]l . Further detailed discussion for SH waves 
is provided in ^ 

More generally, the absence of generating functions for Sm — means that some elements 
of the Willis coupling tensor Sijk should be set to zero in order to complete its definition. The 
relevant elements are necessary to determine stress and momentum but do not enter into the 
equation of motion (j3ip and the sextic system ([2]) with (j2ip because the purely imaginary property 
of the coupling tensor means that mi{Siji — Snj) are the elements of Sm + S+ . Consider the two 
dimensional situation with indices taking only two values so that, on account of the symmetry 
Sijk = Sjik-, there are at most six independent elements. Four of these may be found from (j23p i . 
and one more follows from (126^ 9 using the symmetry property. The single element m * oiY^m IS 
undefined and may be set equal to zero. In the three dimensional situation all but four combinations 
of the 18 independent elements of Sijk are obtainable. Taking an orthonormal triad {mi, n, 1112}, 
the following elements of the coupling tensor are not defined by the effective medium equations 
and are therefore set to zero: iiiq • Smc™/3 + ™/3 ' SmQ™cn a,(3 £ {1,2}. To be explicit, let n lie 
in the 2-direction, then the Sn equation ()23P i defines the nine elements S2jk', these combined with 
the Sm + Sm equations ([26l) 9 yield 5*112, 5*132, 532, and the Sm + S^ equations with the above 
prescriptions give 5ii3 = -5*131, ^313 = -5*331, 5*111 = 0, ^333 = 0. 

3.3.2 Expansion of the Willis parameters 

Explicit insight into the structure of the Willis parameters can be gained from their expansion 
obtained via the Magnus series for the wave number matrix K. In view of (jl4p and (|16p . Ki (0) 
is imaginary and expands as Ki (0) = X^mK^*"^ (0) with odd m and K^™^ (0) ~ iuj"^'^^, while 

K2 (0) is real and expands as K2 (0) = (N2) + Y.m ^2^^ (0) ^i^li even m and K^"^ (0) ~ w™. This 
confirms that Sn is imaginary and vanishes at w = 0. It is easy to check that the right-hand sides of 
(j23k ) and (plj) are zero at kx = and at u, kx = 0, respectively. Based on the forms of Kj (w^, kx) 
as generated by the Magnus expansion, and evident from p4p for the leading order contributions, 
it may be demonstrated that ()26p 9 is consistent with imaginary Sm that is zero at w = 0. It is 
noteworthy that keeping the density p^"^^^ as {p) or as any other scalar is generally not possible 
since this would contradict the pure dynamic term on the right-hand side of (I24p . Finally, it is 
emphasized that, by virtue of (jl7p . the stress-impulse tensor defined as a pure dynamic quantity 
Sijk = Sijk {^) vanishes in the case of a unit cell with any symmetric heterogeneity profile Q (y) 
(regarding the "inaccessible" part Sm — S+j, see ^3.3.ip . 

Equations (I23p and (I26p with polynomials Kj [ijj'^,kx) given by the Magnus series ()13p imply 

that the elastic moduli cl'flj are rational functions of , kx while the density and coupling terms 
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p(off))^ ujSijk are functions of w^, defined by the series 

m=l,2,... 

with real p'^'"-' and imaginary S^Jj^ proportional to a;*", c|^] real or imaginary depending as m is 
odd or even, respectively. These series are similar to ()13l) in that they are majorised by the power 
series in long-wave parameter e. The Magnus series with M terms enables finding M terms of the 
series (fMl) . It is apparent from (f23l ). ([26li ) and ([3]) that Sn and p^'^^)) depend only upon N2 and 
p, thus 

Sn {io) = -iu: (N2)"^ {a^^^ + u^^'^ - af (Ns)"^ a^^^) + ...}, (35) 
p(^«^) (..) = (,) I -.;^(a(^)-a(^) (N^)"^ a«) + . . . , 



with 



(1) T 
2 2 



■ / /' (pN2 - N2p) (= a(^)^), 
Jo Jo 

{af \ a^')(= a^''^)} = ^ ( f H {2N2PN2 - N2N2P - PN2N2, ppN^ + N^pp - 2pN2p}, (36) 

D Jo Jo Jo 

.3 ' = ^ / / /V (2pN2pN2 - 2N2PN2P + N2N2PP - PPN2N2) , 

D Jo Jo Jo Jo 



in which d^, d?!, ... are suppressed (as kept tacit hereafter) and dependence of co- factors on the 
successive integration variables ^, 51, ... is understood. The remaining part of the coupling tensor 
is only obtainable through the combination Sm + S+ , and it depends upon Ni, N2 and p, with 



S„, + S+ = iw( (Ni)T (N2)-^ aW - aW (N2)-^ (Ni)^ - a^) + 0(a;3), where (37) 



For example, the expansions of the inertial quantities for the transversely isotropic layered 
medium discussed in §3.3.21 are p^°^^ = (p) I + p*-^-* + • . . 



^22 = '^^T{c^,Y^ I' j\c^^p - pc^^), (38) 

s<;. . /kv - f (d - If). - .(1 - i)) }. 

and 5'2xi = 5'233 have respectively the same form as p^22 ^ ^222 with C22 replaced by cge in 

1,2- 



3.3.3 Effective medium defined from the Floquet dispersion 

Modelling a dispersive effective medium may be based on a more relaxed approach that abandons 
fitting the matrix to the coefficients of sextic system of wave equations and deals instead with 
the asymptotic secular equation for the eigenvalues iK^ or e*^"^ of or M(T, 0), which is a 
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dispersion equation for the onset of fundamental Floquet branches Ka (w, k^) or uja {kx, K) analyzed 
in [SIEII?]. This gives the same secular equation as that for the iK matrices, and hence preserves 
the long-wave Floquet dispersion but not the displacement-traction vector Wq at the period edges 
(see (fTO|) ). By not fitting all of the physical properties, this type of approach to homogenization 
modelling introduces extra degrees of freedom. In particular, a "modified" effective medium may 
be defined that is asymptotically similar to iK but the matrix Qcg has no pure dynamic terms in 
the diagonal blocks, and hence matches the Stroh-like form ([3|), ^ (though now with (j27P i ). i.e. 
satisfies the standard form of the governing equations ([T]) with dispersive effective coefficients. 
For instance, in the ID case k^ = 0, the matrix 



has asymptotically (to the order of this matrix itself) the same secular equation as the matrix 



The latter "skips" (by construction) the Willis coupling tensor and leads to the same definition 
of the matrix of second-order elastic coefficients (nn)^^^ as in l\23\\ . while the second-order density 
matrix p^^^ = io'^{8L^^'^ (N2)~"'^ — a^^^) following from (j40p is generally different from p^^^ in ([35]) due 
to non-commutativity of (N2) and sl^^ ■ See also the SH example in Appendix B. 

4 Effective medium coefficients for SH waves 
4.1 The wave number matrix 

Consider SH waves in an isotropic medium with periodic density p (y) and shear modulus p (y) . 
The SH state vector ry (y) = [A, iF)^ , where A and F are the amplitudes of u = ^3 and CJ23 (the 
indices correspond to u || X3, n || X2, m || Xi), satisfies Eq. ^ with the system matrix 



The 2x2 case leads to some simplifications not available for higher algebraic dimensions. In par- 
ticular, the two eigenvalues of the monodromy matrix M(T, 0), which are the inverse of one other 
(since detM = 1 due to the isotropy), are defined by the single quantity trM(T, 0). The implica- 
tions are explored in [T^ and only the necessary equations are cited here - the reader is referred to 
[15j for details. The main result is that the wave number matrix, and hence the effective system 
matrix Qefj (w) = iK has semi-explicit form, 




(39) 




(40) 




(41) 



Qoff = i 



( 



K3 -Ki 



) 



sin KT 



K 



[M(T, 0) - I cos KT] , KT = cos 



-1 



(-trM(r,0)), (42) 



where Re cos € [0, vr], Imcos 



> 0, and ±i^r (no subscript) are the eigenvalues of K. 
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4.2 Willis equations and effective coefficients 

Following the general formalism of J}3]tlie effective material is assumed to have constitutive equations 
described by the Willis model, which in this case has only a single momentum component and 
the usual stress components for SH waves in elasticity. Noting that = 0, on account of the 
transversely isotropic axis n, we have 

/ (eff) (cff) 

= ctP 4f ) 543 I I - 2 1 . (43) 

V S43 




These constitutive relations imply, using (jl9li ) , that the governing equation for the SH displacement 
is of the form 

4fn" + (.V=*^)-A:^4f)n = 0> (44) 
where ' means d/dy. The coupling term S43 is absent from the equation of motion, as expected 
from the Willis equations (jl9p for a scalar problem. At the same time, (j43p leads to the state-vector 
system matrix in the form 



Q 



eff 



— C44 lfcxC45 — UJOis) — C44 

,i.2>ff) , ,2^(cff) , >ff)-Vi. ,,Q ^2 Jeff)-!., (eff) ^ ^ 



(45) 



where c'^^^^ = c'^^^ = —cf^^ has been used. 

Setting Qeff of the Willis model equal to that of ()32|) gives the material parameters 

ctP = k~^Ks + K^'K!+u;'p(^^)), c^^ = k~'{K^'K, - K^Ho)K,{0)), (46) 

with Kj{0) = Kj{u}, 0). These may be expressed directly in terms of the elements of the monodromy 
matrix, using the form (j42p along with detK = —K^, 

(eff) Sini^r K^jO) (eff) M4(0)-Mi(0) 

cfP = fe- {.^p^^-^ - K^tP) , = (.^43 + ^^^) , (47) 

where Mj = Mj{T,0) are functions of uj and k^, and (0) means evaluated at k^ = 0. Note that 
the expressions for p^*^^^ and also follow from the equation of motion (j44p and its solution 
u{y) = u{0)e'^y, using k^ = for p^^^^uj). 

Explicit expressions for the low-frequency long-wave expansion of the material parameters may 
be found in the same manner as in §3.3.21 for the general case. The starting point is the Magnus 
expansion 

Qeff = (Q) + iK(i) + iK(2) for the SH wave number matrix. Details of the analysis and 
a summary of the results are presented in Appendix B. 



4.3 Examples and discussion 
4.3.1 A bilayered unit cell 

The general formulation is illustrated by the case of a two-component piecewise constant unit 
cell. Specifically, consider a periodic structure of homogeneous isotropic layers j = 1,2, each 



12 



with constant density pj , shear modulus pj and thickness dj . The monodromy matrix M (T, 0) 
gQ2rf2gQirfi = ]y[ (^^^ j^^-j -j^Qg ^jjg well-known form 



cos ■(/'2 cos V'l - ^ sin ^^2 sin V'l cos ?/'2 sin^i - sin ^^2 cos V'l \ , . 

-Z7i cos'(/'2 sin'(/'i — Z72 sin'(/'2 cos^/'i cos '(/'2 cos V'l — ^ sin -02 sin i' 



where tpj = djJp.^ ^Pj^"^ ~ is the phase shift over a layer and 7^ = fijipj/dj, see [15]. Figures [T] 




Figure 1: The effective material properties of the bilayered SH case for kx = 0: elastic moduli, 
inertial parameters and the effective wave number are plotted in blue, black and red, respectively. 
The frequency range includes the first band edge which is at the frequency where JieK = vr first 
occurs. Only the real parts of the quantities indicated are plotted. For frequencies in the stop band 
the imaginary parts of 044^^ p^^^"^ and K are non-zero but not shown. 

and [2] show the computed parameters for the case of layers of equal thickness, di = d2 = 1/2, with 
pi = 1, ci = 1; p2 = 2, C2 = 2, where Cj is the shear wave speed (c^ = p/p). Figure [T] shows the 
effective parameters for propagation normal to the layers {kx = 0). The vanishing of both c^^^^ and 
p^^^^ at the band edge at w = wi f» 2.6 is expected on the basis of the fact that Qefr is singular at 
the band edge and scales as {uj — uji)~^/'^ near it [15]. Referring to the 12-element in Eq. (j45p . this 
implies first that c^^^^ oc (w — wi)^/^ and then, from the 21-element and the finite value of det Qefr, 
that oc - The square root decay of both clf ^ and p'^^^) is apparent in Figured! 

The wavenumber is finite, /c^; = 1, in Figure [2j This has the effect of increasing the frequency 
of the band edge, and introducing a range of frequency from u; = up to the cut-on at w ~ 1.7 in 
which the effective wave is non-propagating. Note that p^^^^ and S43 are unchanged from Figure 
[1] while the elastic modulus c~l^'^ is different, and tends to zero at the new band edge as expected. 
The non-zero kx leads to non-zero c'l^\ and the parameter c~^^^ becomes complex- valued at the 
kx = Q band edge. Only the real parts of the quantities are shown in both figures. No attempt 
is made here to discuss their imaginary components, which requires careful analysis of the branch 
cuts and is a topic for separate study. 
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Figure 2: The same as in Figure □ but for = 1.0. The additional parameters 4f ^ and cjf ^ are 
relevant to 7^ 0. Only the real parts of the quantities are plotted. 

4.3.2 Reflection and transmission of a half-space of effective material 

As an example of the type of boundary problem that can be solved using the effective medium 
equations, consider reflection-transmission of SH waves at a bonded interface y = between the 
half-space of the periodically stratified medium (y > 0) and a uniform half-space {y < 0) of isotropic 
material with po, fiQ and cq = \/ po/po- A SH plane wave is incident from the uniform half-space 
with propagation direction at angle 6 from the interface normal. The total solution is taken as 

u{x,y) = e'^^"" X <^ L / }, U - , ^. , ^ ^ -(sin0, cos^). (49) 

The reflection and transmission coefficients R and T follow from the continuity conditions for 
particle velocity and traction at the interface. They may be expressed in standard form using 
SH impedances defined as Z± = — o"23/u|y=o± ■ The impedance in the uniform half-space is Z_ = 
PqCq cos 9. The impedance for the effective medium follows from (143p as 

-1/ 



cif + fc.4f -a;543). (50) 



Z+=uj-^{Kc 

This is identical to the impedance of the periodically stratified half-space because they both imply 
a ratio of components of the outgoing eigenvector w which is common to M(T, 0) and K. In these 
terms, the continuity conditions for displacement and traction yield the exact result 

Z — Z 2Z 
1 + R = T, Z_{l-R) = Z+T, => E=___±, T=^-^. (51) 

Figure [3] shows \R (uj) \ and \T (u) \ calculated for normal incidence from a uniform half-space with 
pi = 1, ci = 1 on a periodic structure of two layers with = 1, ci = 1, /32 = 2, C2 = 2, which was 
used in Figs. 1 and 2. As expected, \R\ < 1 with total reflection in the stopband. 
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The explicit dependence of the reflection coefficient on the effective medium parameters c'^^^ 



-45 



and 5*43 means that, in principle, measurement of R via experiment can provide useful knowl- 



edge for their determination. 
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Figure 3: The magnitude of the reflection and transmission coefficients of (I5ip for normal incidence 
(^ = 0). 



4.3.3 Uniform normal impedance 



It is instructive to consider the particular case oi kx = Q with z = y/p{y) fi (y) = p{y) c (y) 
independent of y, i.e. z = zq . The 2x2 matrix Q (y) is then a scalar multiple of a constant matrix, 
and Mi(0) = Mi{0), M2(0) = (iujzo)-^ sinKT (where (0) stands for = 0). As a resuh, by 
(f46l) . the effective parameters at any w retain their values obtained from static averaging: 543 = 0, 



C44 (a;,0) = {p~^)~^ and p^'^^^Lo) = {p). This simplification is a consequence of the fact that 
constant z implies constant eigenvectors of the SH matricant M (y, 0) and hence no refiection of SH 
waves normally propagating through a periodic structure, which is in accordance with the physical 
meaning of the impedance z. Consistency is also observed in that the effective impedance z^^^^ 
defined through the above effective parameters is equal to zq, 



Zq. 



(52) 



The only effect of the inhomogeneity is to speed up or retard advancing waves according to the 
effective speed c^^^^ = z^^^^ / p^'^^\ which in the present case follows from (I52p as c^®^^ = zq/p^'^^^ = 



4.3.4 Discussion 

In the case of purely uni-dimensional motion, kx = 0, the system (|43|) involves only the first three 
parameters of Eq. (j37|) : c^^^ , p^*^^) , and S'43 . Willis [H] derived expressions for the same quantities 
for a laminated medium. His coefficients |14l Eq. (3.30)] relate weighted means of strain and 
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velocity {{we), {wu)) to ensemble means of stress and momentum density {{cr), {p)), where w is a 
general weighting function first introduced in [13], such that ensemble means correspond to w; = 1. 
The Willis parameters derived here, e.g. ()43p . concern strain and velocity at the single point y = 
in the unit period, and therefore correspond to the specific weight function w{x) = 2L6{x) in the 
notation of [14] . It is important to note, however, that the stress and momentum density used here 
are not ensemble averages but are quantities associated with the same point in the unit period. 
This identification, for instance, means that the solution of the reflection-transmission problem of 
§4.3.2l is in fact the deterministic solution. In summary, while the governing equations are the same 
in both cases, the Willis parameters developed here do not bear a one-to-one correspondence with 
those in [T4] . 

5 Conclusion 

A fully dynamic homogenization scheme has been developed for periodically layered anisotropic 
elastic solids. In the process, the dispersive and nonlocal Willis model has been shown to provide 
an optimal constitutive setting for the effective medium. The crucial point of the present method is 
the insistence that the matrix of coefficients Qog of the sextic system of elastodynamics equations, 
for whatever homogeneous effective medium is considered, must exactly match the Floquet wave 
number matrix K of the periodic system. This is not a low-frequency long-wave approach, so long 
as K is defined at the given frequency cj and horizontal wave number k^- The wave number matrix 
K = K{uj,k x) is an analytic function of lj and kx that may be explicitly defined via the Magnus 
series expansion, which is guaranteed to converge below the first Floquet stopband at the edge of 
the Brillouin zone. The choice of constitutive model for the effective medium is critical. We have 
demonstrated that the standard anisotropic elasticity theory does not suffice as it cannot provide a 
Qeff to properly account for dynamic terms appearing in the wave number matrix K {to, kx)- On the 
other hand, the Willis model for the effective medium, which includes coupling effects, can allow 
us to associate elements of the effective system matrix Qeg with elements in K. The main results 
are contained in Eqs. (j23p and (I26p which infer the material parameters of the effective Willis 
medium from K. Invoking the Magnus series, explicit expressions for the low-frequency long-wave 
expansion of these effective Willis parameters have been found and the accuracy for their truncated 
asymptotics has been estimated. 

The example of SH plane wave reflection and transmission considered in ^4.3.21 indicates the 
type of application possible using the Willis effective medium. The point is not so much to provide 
new solutions for layered media, although it is simpler to formulate and solve such problems using 
equations for a homogeneous model. The potential power of the dynamic effective medium model 
is that it is possible to relate the effective properties of the Willis material to measurable dynamic 
quantities. Thus, the reflection and transmission problem illustrates how the reflection coefficient 
R depends on a certain combination of the Willis parameters. Measurements of i? = R{uj,kx) 
provide a means to characterize periodic layered systems as equivalent homogeneous but dispersive 
materials. Other problems that may be considered are, for instance, surface wave propagation in a 
periodically layered half space, waveguides comprised of periodic layers, and point forces. 
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Appendix 

A. Estimates for the long-wavelength expansion of K 
A.l Auxiliary notations 

For any 2d x 2d matrix A consisting of d x d blocks Aj (J = 1, ...,4), denote a 2x2 matrix of 
matrix norms ||Aj|| by ||A||2x2 : 




Note that ||AB||2x2 < ||A||2x2||B||2x2i where and hereafter a matrix inequality is understood as 
that between the corresponding matrix elements. Let 

= {Nif + {N2) {N3) , q{y)= max (||Nj (y) \\/{Nj), p {y)/{p)) , (54) 

where (•) is the averaging symbol and (Nj) is an average of a norm Nj (y) = ||Nj (y) || of the 
Stroh-matrix block (J = 1,2,3). (||Ni|| = IIN4II). Both u and q{y) are physically dimensionless 
and strictly positive. The magnitude of is usually of the order of 1 unless a high-contrast case 
affecting the averaged profile; in turn, the averaged value (q) , strictly speaking, satisfies 1 < (g) < 4 
but is typically close to 1 as well. Let us define the long-wave small parameter e as 

e = kT withfe = V^2+^2/y2^ V^={{Nif + {N3){N2)){{p){N2)y\ (55) 

It is noted that e = z/e may equally be taken as such a parameter and that it can actually be 
replaced everywhere below by a smaller value 

? = Tmax (fc, (iVi) , k,,^{N2) (A^a), ^V(p) {N2)) < ve. (56) 

Next, define 

C - („ imi) . -3 (.) . CQ to) ^ . ^^^^ _ J , ,5T) 

where Q (y) is given in ([3]), and the normalization of Q (y) provides a common estimate ||Qj (y) \F < 
ueq (y) for the blocks of Q, from which follow the estimates for ||Q,/ (y) ||. Denoting the (constant) 
2x2 matrices 

H = (I I) > I|C||2X2 - C = (J _^^°^^^) , n = {q) C-HC, (58) 

and using the notation (|53p enables us to write the blockwise estimates for Q (y) and for Q (y) in 
the form 

l|Q(?/)|l2x2^<-^9(2/)H ^ ||Q(y)||,,, = C-iQ(y)||,^,C<^n 

^I|(Q)|I2X2<(||Q|I2X2)<"- 

Note the meaning of as a matrix of upper bounds of blocks of the statically averaged system 
matrix (Q). 



17 



A. 2 Estimates for the Magnus series iK = Ylm=o^^^"^^ 

An elegant proof that the Magnus series (jl3|) converges for ( 1 1 Q 1 1 2 ) ^ < IE] is somewhat imphcit 
in that it does not provide fully explicit estimates of the series terms and remainder. These may 
be obtained for a narrower range by adapting the derivation detailed and referenced in It 
proceeds from the estimate 

||K("^)||r<7r(er(||Q||))™+\ e = 2( r[2 + x(l -cot x)]-idx)"^ = 1.8400... (60) 

Jo 

with specifically the matrix norm || • II2 as kept tacit below. (Note aside that the physical dimension 
of (QT)" is the same for any n as that of QT and K'^"*)T.) Let iK = X]m=o ^K^™) be the same 
series (fT3]) but related to the matrix Q in place of Q. Applying (i60]l yields the blockwise estimates 
for series terms K^"^) and, hence, K^*") in the form 

||K^"^||2x2T < 7rX"^+^H ^ ||K(™)||2x2 < where X = ^ue {q) . (61) 

Assume hereafter that A < 1 (which is within the convergence radius A < that follows from 
the result (||Q||)r < vr of [S]). By dSU, the residual series R(^^) = Em=M^K("") and R(*^) = 
E™=M^K(™) satisfy 

||R^^^^||2x2r < ^(1 - X)-^X^'^+^H ||R(^^ 112x2 < vre(l - A)-iA*^n, (62) 

so that they decrease as M grows and tend to zero as M — )■ 00. 

Knowing the upper bound of HR^^^-* ||2x2 evaluates the sufficient number of terms to be kept in 
the Magnus series to ensure a desired accuracy of truncation for a fixed long- wave parameter e, or 
else provides the value of e that ensures this accuracy for a given truncation step. The accuracy is 
gauged by the matrix fl of blockwise bounds of (Q) , see (f59]) . For example, let the Magnus series 
(jl3p be truncated as iK = (Q) + iK^^^ + iK^^^ and the remainder R^'^^ discarded. According to 
(f62l) o. ||R(^)||2x2 < 5.8A3/(1- A). Thus taking the spectral range as ue {q) < 0.24 or < 0.128 
ensures ||R^^^||2x2 < or < O.lft, respectively (note that truncating zK by (Q) at i^e (q) < 0.128 
discards ||R(^)||2x2 < 1.79^). 

It is noted that the above mentioned sufficient criterion for the Magnus series convergence 
and the bounds for its terms restrict uj and kx by imposing conditions on the norm ||Q|| of 
Q = Q (y; '^i ^x)- At the same time, if Q is independent of y then the Magnus series certainly 
converges at any u and kx- Hence if the inhomogeneity is relatively weak so that ||Q— (Q)|| is 
markedly smaller than ||Q|| , then the uj and kx bounds on the Magnus series can be eased by using 
a different approach that is based, instead of ([60|) . on the estimate IIK^"*)]] < Cm{\\Q— (Q)!!)'" 
(C > is some constant). The latter acquires a growing factor m but is hinged explicitly on 
HQ— (Q)||, which is why this approach can be much more advantageous for small ||Q— (Q)||- 

A. 3 Invertibility of K2 

Derivation of the effective material constants requires inverting the block K2 (see §3.2). Evidently 
iK2 = i (N2) + R2^'' is assuredly negative definite (like (N2)) if || (N2)~"^ R2"^'' || < 1. Inserting 
||R2^''|| ^ TT^ {q) (A2) A/ (1 — A) from (j62|) 9 and resolving the resulting inequality with respect to 
ue {q) yields 

ue{q)<{2 + n{q){N2)\\{N2r'\\)~\ (63) 
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This condition can be further improved by using more precise estimates for the low-order terms of 
the Magnus series. For this purpose, it is suitable to proceed from K^™"^ defined by (I13p with Q. 
Using f^' ... /|-- q (?) ...q {^^) = ^^^^H and HK^'")!! = ^HK^'")!! gives for m = 1, 2 : 

(64) 



(N2)"^R^^^|| < 1 if X < min (0.41, (l + 1.4(g) (iVs) || {^2^^^ 



||K«||<z.e((7)(iV2), \\K^?\\<lue{qf{N2), 
where the right-hand sides are smaller than in (j6ip with m = 1,2. Exploiting (|64p leads to 

II (N2)-^r(^)|| = II (N2)'^ (eLi^k(-) +R(3)) II 
< 2 {q) {N2) II (Ns)"^ II (0.272X + 0.132^2 + 2.93XV(1 - X)) (65) 

^2) 

Hence, by ([65]) . taking a typical value {q) {N2) \\ (N2) II ^ 1 guarantees that K2 is invertible in 
the spectral range X < 0.41, i.e. ue {q) < 0.22. 

A. 4 Estimate for the eigenvalues of K 

The upper bound for eigenvalues of K may be evaluated via ||i^Q|| < ||K||; however, taking 
note that are also the eigenvalues of K yields a better estimate with regard for ()62p 1 as follows 

ll-f^all < ||K|| = ||r(°)|| < 2||R(°)||2x2 < 27rr-^(l - X)~^X (66) 

(note that here may certainly take either real or complex values). For instance, ||ira||r < 0.62 
in the long-wave range i^e {q) < 0.128 {X < 0.236) that was shown in §A.2 to ensure accuracy of 
truncation of K after K^^^. Note that vanishing of the right-hand side of (|66p at X = 0, i.e. at 

e = Ty^k^ToJyV^ = (see ([55])) is in agreement with ^ = (^^ ^^^^^ at a; = 0, = 0. 
A. 5 Convergence of dK./dkx 

Consider {dK/dk^)i^^^Q = K' (0) as defined by Eq. ([29]). The series ([MIl) and ([M) assuredly 
converge if so does the same series for InCMC"^ = C (InM) where C is introduced in (fSTli V 
In turn the series for InCMC"^ assuredly converges if 

||CMC~i-l|| < ||exp(||Q||)-l|| < exp(2z/e(g)) - 1 < 1 ^ 2iye{q)<ln2, (67) 

where it was used that ||Q|| = 2ue (q) due to dM) and ||H|| = 2. For M = M (0) with k^ = 0, 
which is the case in hand, the condition (f67|l reduces, using (|56|) . to 



2?(0) (g(0)) < ln2, (68) 



where ?(0) = rmax(^^(p) {N2)), g (0) = max (||N2 (y) ||/(iV2), p (y)/(p)) {1 < {q (0)) < 2) . 

The integrals in (j28b ) and ()29b ) exist provided the matrix I -1- x(M(0) — I) is invertible for any 
X € (0, 1). According to [16], this is guaranteed if 

2?(0)((?(0)) <7r, (69) 

which implies that the eigenvalues e^^°'^^^'^ of M (0) do not attain real negative values. This 
condition coincides with the sufficient condition for the convergence of Mag nus scries cit hx — 0. 
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B. Effective medium coefficients for SH waves at long wavelength 

According to (fT3]) truncated by K^^^ , 



Qcff 



K 



(1) 



if,) - {p) co' + iq 



(2) 



(1) 



(70) 



where ^^2^3 = 0, KY'' = and (omitting the integration variables) 



.(2) 



iT 



1 /•? 



"'0 



On the other hand, Eq. (j45p taken to the same order as (j70p and written with Voigt index notation, 
yields 



Q 



where c 



(1)+ 
45 



eff 



-54 



{p ^> (fcxcis 



45 '^'-'43 



+ 



{p ^) (^k^c 



(1) 

45 "^'-'43 



c\J has been used, and in addition c, 



JO) 



0, pW = with, as noted earlier, 5^3^ = 0. From (f70]l -f72 



■44 



' '-55 



ik-), ^4 



(72) 



.(1) 
-55 



^45 



.(2) _ 



-44 - \k V ' ^^^5^ - p'''^'' = Kr - {p~t ^ > 

which identifies 045^(^3;) = —c\^^{kx), c'l} (uj'^,k'^^ and 543^(0;). The density correction p'^'^\< 

follows from the final identity as equal to the expression (f38|) i with C33 — )■ p, and c"^^ (a;^,fc^) is 
then uniquely defined. Similar results follow for the SH waves in a monoclinic periodic medium. 

Note that the secular equation of the SH-wave matrix Qcg given by ([70]) can equally be associ- 
ated, up to the same order as Qefr itself, with other matrices such as, e.g., the matrix 



Qeff 







Sp)kl-{p)u^ + Kf^-{p~^y'K[ 



R-(l)2 



(74) 



.(1) 
-45 



which has zero diagonal components like in (j4ip and thus leads to Eqs. (j43p . (j73p but with zero c 
and 543^. However, as remarked in ^3.3.31 the effective polarization of the SH displacement-traction 
eigenmodes is asymptotically defined by the eigenvectors wi^2 of (truncated) Qcfi = ^K, while its 
definition from the eigenvectors wi^2 of Qcfr is different since wi^2 7^ "^1,2 to the first order in e. 
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